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1 We present a statistical representation of the aggregate effects of deep convection on the 

2 chemistry and dynamics of the Upper Troposphere (UT) based on direct aircraft 

3 observations of the chemical composition of the UT over the Eastern United States and 

4 Canada during summer. These measurements provide new and unique observational 

5 constraints on the chemistry occurring downwind of convection and the rate at which air 

6 in the UT is recycled, previously only the province of model analyses. These results 

7 provide quantitative measures that can be used to evaluate global climate and chemistry 

8 models. 
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1 Deep convection is a highly efficient mechanism for the vertical transport of air from 

2 near the Earth’s surface (0-2 km) to the Upper Troposphere (UT) (6-12 km) (1-5). 

3 Typical convective storms have spatial scales of tens of kilometers and vertical velocities 

4 as large as 15 m sec ' 1 (6), making their local influence in the UT extremely strong. The 

5 rapid upward flow is balanced by much slower descending flow that occurs over a larger 

6 spatial scale (7). Convection is also associated with lightning a significant source of NO x 

7 (NO x = NO + NO 2 ) in the UT (8, 9). The source strength and spatial distribution of 

8 lightning NO x emissions are not well known, with estimates ranging from 2-20 Tg(N) yr ' 1 

9 for the global average (10), compared to 25 Tg(N) yr ' 1 from fossil fuel combustion (11). 

10 Although there have been a number of case studies of the chemical effects of individual 

1 1 storms (12), studies of the aggregate effects of convection on the chemical composition 

12 and radiative forcing of the UT have been largely the province of modeling and theory 

13 (13, 14). Here we describe measurements that provide a direct link between an 

14 observable property and the ensemble of convective events. 

15 

16 The chemical and radiative consequences of convection and lightning are known to be 

17 large (2, 3, 15). Upper tropospheric O 3 , either transported directly from the boundary 

18 layer via convection or formed in situ following detrainment of convectively lofted O 3 

19 precursors (NO x , odd hydrogen radicals (HO x ) and hydrocarbons) in the outflow region, 

20 directly impacts climate through a positive radiative forcing (15). Additionally, deep 

21 convection accounts for a significant fraction of the net flux of moisture from near the 

22 Earth’s surface to the UT (16). Thus, the rate at which the UT is turned over by 
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1 convection has important implications on the hydrological cycle and the magnitude of the 

2 water vapor feedback on global temperature (7 7). 

3 

4 In this study we describe a method for calculating the time air spends in the upper 

5 troposphere following convection from in situ measurements of the chemical composition 

6 of the UT and discuss the chemistry occurring in the outflow region as a function of time 

7 since convection. We use measurements of NO 2 ( 18, 19) (NO x is calculated from NO 2 , 

8 O 3 , HO 2 and photolysis rates), HNO 3 (20, 21), OH and HO 2 (22), O 3 (23), aerosol number 

9 density (24), actinic flux (25), CO (26) and CO 2 (27) obtained during the Intercontinental 

10 Chemical Transport Experiment - North America (INTEX-NA) aboard the NASA DC -8 

1 1 (28). Measurements were made at altitudes between the surface and 12.5 km, over a wide 

12 area of the US and Canada, west of 40° W and between 30 and 50° N. There were a large 

13 number of vertical profiles allowing a reasonably unbiased statistical sampling of air over 

14 this region during July and August of 2004. 

15 

16 We use the deviation of the observed NO x to HNO3 ratio from steady-state as an indicator 

17 of convective influence. The NO x to HNO 3 ratio is reset to near infinity in moist 

18 convection as a result of preferential wet scavenging of HNO3 relative to NO x (i.e., the 

19 Henry’s Law Constant for HNO 3 is ~ 10 8 times larger than for NO x ) (29). Further, 

20 lightning initiated NO x production, often coincident with convection, dramatically 

21 enhances NO x in the outflow region. The coupling of these processes makes the NO x to 

22 HNO 3 ratio in the UT an effective indicator of convective influence, where NO x /HN 03 

23 » 1 is indicative of recent cloud outflow (30, 31). In the days following convection, the 
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1 ratio decays toward steady-state providing a chemical clock that marks the time an air- 

2 mass has spent in the UT following convection (32). A number of previous studies have 

3 used species which have no upper tropospheric source (e.g. CH 3 I) (33), or alternative 

4 chemical ratios to provide estimates of age of air in the UT {34, 35). Our study is unique 

5 because of the availability of high time resolution NO 2 , OH and HNO 3 measurements that 

6 allow us to build a much more extensive data set than previous studies. After the initial 

7 turbulent mixing in the near field of the convection, mixing is slow, thus the time 

8 evolution of N 0 X /HN 03 following convection depends largely on the partitioning of NO x 

9 (between NO and NO 2 ), the concentration of OH and the actinic flux. 

10 

1 1 Reactive Nitrogen Partitioning in the UT 

12 The only significant chemical sinks of UT NO x are reaction with OH to produce HNO 3 

13 (Equation 1) and nighttime loss through NO 3 (Equation 2a-b) followed by hydrolysis of 

14 N 2 O 5 to produce HNO 3 {36). NO x is regenerated by nitric acid photolysis (and 

15 subsequent NO 3 photolysis to NO 2 ) and reaction of OH with HNO 3 (Equations 3 and 4). 


16 

N0 2 + OH HNO, 

( 1 ) 

17 

NO 2 + 0 3 — ) 

>no 3 +o 2 

( 2 a) 

18 

NO, + NO 2 

-> n 2 o 5 

( 2 b) 

19 

HNO- s hv 

>oh+no 2 

(3) 

20 

hno 3 + oh ^ 

>no 3 +h 2 o 

(4) 


21 Including the altitude dependent rain-out rate for HNO 3 (k ra i n -out) as derived by Giorgi and 

22 Chameides (37), the expected steady-state N 0 X /HN 03 is: 

23 
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2 

3 Our observations show the NO x to HNO3 ratio to be much higher than the ratio described 

4 by equation 3 at altitudes greater than 6 km (c.f. Fig. 2). The difference between the 

5 observed ratio and that predicted by equation 3 grows with altitude, reaching a maximum 

6 at 10 km. Previous observations of NO x and HNO3 (either measured directly or 

7 calculated from observations of NO x , PAN and NO y ) have shown the N0 X /HN03 ratio to 

8 be significantly larger than the steady-state prediction in the UT (30, 31, 38-42). This has 

9 been shown to be primarily a result of convection and lightning reinitializing the system 

10 before steady-state is achieved (30, 31). Although a series of other hypotheses have been 

1 1 put forth (40-42), we (like Jaegle et al.) find no evidence for a mechanism other than 

12 convection responsible for holding NO x /HN03 out of steady-state in the UT. 

13 

14 Chemical Signatures of Convection 

15 Fig. 3 depicts one of many convectively influenced air-masses sampled in the UT during 

16 INTEX-NA. Three distinct convective events (40-80 km wide) are identified by 

17 enhancements in NO x /HN03 in Fig. 3a. Coincident enhancements are present in SO2, an 

18 indicator of a recent boundary layer source for this air, and Ultra- fine Cloud 

19 Condensation Nuclei (UCN) (3 < D p < 10 nm), an indicator of cloud detrainment (Fig. 

20 3b) (43, 44). Sharp decreases in CO2 also indicate the convective lofting of boundary 

21 layer air depleted in CO2, a result of photosynthetic activity (Fig. 3c) (44). 

22 Enhancements in CO, CH2O and various hydrocarbons, relative to the surrounding UT 
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1 air, were also observed in these plumes, further indicating that these parcels originate 

2 from the Planetary Boundary Layer (PBL). Backward air trajectories, initialized along 

3 the flight track, coupled to the spatial and temporal distribution of cloud-to-ground 

4 lightning strikes, indicate that the sampled air-mass was recently influenced by lightning 

5 approximately one day prior to DC-8 sampling (c.f. Fig. 3 bottom panel) ( 45 ). Such 

6 features with high N0 X /HN03 were observed throughout the UT during INTEX-NA. 

7 

8 To assess the extent to which the UT over the Eastern U.S. and Canada during the 

9 summer of 2004 was influenced by convection and describe the chemical evolution of 

10 convective outflow, we use a constrained time-dependent photo-chemical box model to 

1 1 map the observed NO x /HN03 to the time since the ratio was last reinitialized. The model 

12 is described in the supplemental information included with this article. It is initialized 

13 with observations at 1km vertical intervals from 6 to 12 km. The derived timing indicator 

14 for the convectively influenced air sampled on 1 1 August 2006 is shown in Fig. 3d. The 

15 properties of the ensemble of our measurements are shown in Figs. 4-6. 

16 

17 The aerosol size distribution provides an independent indicator of air recently detrained 

18 from clouds. Cloud-processed air is depleted of aerosol surface area permitting new 

19 particle formation in the outflow region ( 43 , 44 ). Fig. 4a depicts the fraction of 

20 condensation nuclei found in the 3-10 nm bin as a function of time since convective 

21 influence. The fraction of particles in this ultra-fine mode is largest during the first few 

22 days confirming that the NO x to HNO3 ratio, and the timing indicator derived from it, is 

23 reinitialized in the UT by cloud processing. Strong enhancements in CH3OOH/H2O2 (not 
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1 shown), also an indicator of recent cloud processing ( 33 , 46 ), were observed during the 

2 first two days after cloud processing. 

3 

4 As expected, both elevated NO x and suppressed HNO3 are observed at short times (c.f. 

5 Figs. 4b-c). Enhancements in NO x during the first few days is indicative of convection of 

6 boundary layer and/or lightning NO x ( 47 , 48 ). The suppression of HNO3 at short times is 

7 clear indication of HNO3 scavenging during convection. Fig. 4d confirms that reactive 

8 nitrogen (NO y = NO x + total peroxy nitrates + HNO3) is conserved during the chemical 

9 processing following convection, a fact which provides further support for the use of 

10 N0 x /HN 03 as a marker representing time since convection. 

11 

12 Chemical Processing in Convective Outflow 

13 Mapping the ensemble of observations made throughout the UT onto the coordinate of 

14 time since convection allows us to assess the chemical and dynamical processes 

15 occurring following convection, without attempting a Lagrangian convection study. In 

16 this analysis we concentrate on the time evolution of CO and O3. 

17 

1 8 The time evolution of CO, following detrainment into the UT, is set by the abundance of 

19 OH and the rate at which the convective plume entrains air from the background UT (c.f. 

20 Fig. 5a). Due to the direct dependence of the chemical clock on HO x , we constrained 

21 both OH and HO2 to the observations as a function of NO x and pressure in the time- 

22 dependent model used to generate time. As a result, we can iterate the model to 

23 determine the proper mixing rate of the convective plume by matching the modeled and 
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1 observed time evolution of CO following convection. Using this approach for a series of 

2 long lived species (e.g. CO, CH 4 , CH 3 OH and others), we calculate an average mixing 

3 rate of 0.05 ± 0.02 day ' 1 following detrainment into the UT. This is in good agreement 

4 with the upper limit of 0.06 - 0.1 day ' 1 determined by Ray et al. from observations of 

5 convective plumes observed in the stratosphere during the CRYSTAL-FACE mission 

6 (49), however it is slower than the 2 day dilution time-scale determined by Wang et al. 

7 from observations during the SONEX Experiment (32). Since the DC -8 did not routinely 

8 sample in the turbulent environment directly surrounding convective outflow, this mixing 

9 rate likely reflects diffusive and shear induced mixing subsequent to the initial turbulent 

10 mixing occurring during detrainment from the convective system. Entrainment of UT air 

1 1 either during convective lofting or cloud detrainment is discussed in the next section. 

12 

13 The O 3 mixing ratio as a function of time since convection is shown in Fig. 5b. We find 

14 that on average, convectively lofted air-masses contain less O 3 than the background UT. 

15 This result is consistent with the observed vertical gradient in O 3 observed over the 

16 continental US during INTEX-NA, with lower O 3 in the PBL than above (50). Rapid 

17 changes in the O 3 mixing ratio are observed during the first two days following 

18 detrainment, with the observed O 3 15 nmol mol ' 1 above the initial value by the end of day 

19 two. The observed rate of increase slows exponentially with an asymptote at long time 

20 approaching zero and the O 3 mixing ratio approaching a constant value of 85 nmol mol' 1 . 

21 This is a surprising result, as our model of the O 3 rate of change never approaches zero, 

22 but continues to predict a net increase of 3 nmol mol 1 O 3 day ' 1 at the end of day five 

23 (51). 


9 



Bertram et al. Revision 2 30 October 2006 


1 

2 Constraints on the Convective Turnover Rate of the UT 

3 The convective turnover rate of the upper troposphere is critical for accurately describing 

4 NO x , HO x and O 3 chemistry in the UT (52). However, at present there is a paucity of 

5 observation based constraints available (either meteorological or chemical) to test the 

6 aggregate effects of convection in the current generation of global chemistry and climate 

7 models. To determine the convective turnover rate of the UT from the observations 

8 presented here, both the extent to which the UT is influenced by convection and the 

9 fraction of PBL air in the convectively influenced air-masses must be known with high 

10 confidence. 

11 

12 To determine the fraction of PBL air contained in fresh convective outflow, we use 

13 observations of insoluble long-lived species made throughout the INTEX-NA campaign 

14 over the continent. Assuming that we conducted a statistically unbiased sampling of both 

15 the boundary layer and free troposphere during INTEX-NA, we can calculate the fraction 

16 of PBL air present in fresh convection (f) through the following equation: 

17 

1 s LL„.„ = / EL*, + (i - /) LL w 

19 

20 where [X]ut (t=o) is the mean mixing ratio of species X in fresh convective outflow (as 

21 identified using our timing indicator), [X]ut is the mean mixing ratio of species X in the 

22 UT (7.5-11.5 km) and [X] sur f aC e is the mean mixing ratio of species X between 0-1.5 km. 

23 Using observations of CO, CO2, CH3OH, CH4 and C 2 H 6 we calculate the fraction of PBL 
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1 air in fresh convection to be 0.19 ± 0.05, 0.1 1 ± 0.03, 0.26 ± 0.05, 0.15 ± 0.05, and 0.34 ± 

2 0.09, respectively. We calculate a weighted average for the fraction of PBL air in 

3 convective outflow of 0.17 ± 0.02 by weighting each value by the inverse square of its 

4 uncertainty. This implies that convectively lofted PBL air rapidly entrains the 

5 surrounding air either during ascent or in the turbulent environment of the detraining 

6 flow. These results are consistent with: i.) the observations of Ray et al., who determined 

7 the fraction of tropospheric air in convective plumes sampled in the stratosphere to be 

8 between 0.1 and 0.4 (49), ii.) the observations of Cohan et al., who calculate the fraction 

9 of BL air in fresh convection outflow to be between 0.32 and 0.64 from observations of 

10 CHBr 3 , CH 3 OOH and CH 3 I in fresh convection (33) and iii.) the modeling studies of 

1 1 Mullendore et al., who calculate the fraction of PBL air present in the convective outflow 

12 region of a supercell storm to be 0.26, 10 hours after storm initialization (S3). 

13 

14 Fig. 6 a shows the normalized frequency distribution of the observed time since 

15 convection based on the ratio of NO x to HNO 3 . We find that 54% of the air sampled 

16 between 7.5 and 11.5km had been influenced by convection during the past two days. 

17 The convective outflow was strongest between 9.5 and 10.5 km, where the fraction of 

18 sampled air that is less than two days old exceeds 69%. The vertical distribution 

19 presented here is consistent with previous observations and model analyses of convective 

20 outflow to the UT from individual storms (4, 54) and the vertical distribution of 

21 convectively influenced laminae observed in O 3 sonde data from the summer of 2004 

22 over the northeastern U.S. (55). The shift toward longer times between 10.5 and 1 1.5 km 

23 suggests that either convective cloud tops on average do not extend higher than 10.5 km 


11 



Bertram et al. Revision 2 30 October 2006 

1 over the mid-latitude during the summer ( 56 ) or that transport of stratospheric air, rich in 

2 HNO 3 , contributes to keeping the NO x to HNO 3 ratio low at altitudes greater than 10.5 

3 km (57). 

4 

5 To constrain the turnover rate of the UT from the ensemble statistics generated from our 

6 calculated time since convective influence (Fig. 6 a), we constructed a two dimensional 

7 model of the UT. We assume that it takes 4 days for any individual box to pass through 

8 the sampling region and that each box has not been influenced by convection upon 

9 entering the sampling window. Every six hours we: i.) represent convection by 

10 randomly reinitializing the age of x% of the boxes in the sampling domain to 0 (the value 

11 of x is determined by the turnover rate (varied between 0.05 and 0.2 day" 1 ) and the 

12 fraction of PBL air contained in fresh convection (assumed to be a constant at 0.17)) and 

13 ii.) dilute each box with the mean value of the adjacent 8 boxes at the rate of 0.05 day' 1 . 

14 

15 Fig. 6 b depicts the observed and three calculated normalized frequency distribution of 

16 time since convective influence between 7.5 and 11.5 km. The shape of the distribution 

17 suggests that UT air sampled during INTEX-NA was largely influenced by convection, 

1 8 and that convective ly lofted plumes did not have sufficient time to either mix or age prior 

19 to sampling by the DC- 8 , but instead were transported to the East out of the domain. 

20 Frequency distributions of time since convection, using three different convective 

2 1 turnover rates, in the Eastern half of the 2-D UT model analysis (where we sampled most 

22 frequently) are also shown in Fig. 6 b. Assuming the DC -8 made a statistically unbiased 

23 sampling of the continental UT during summer, the best match between the model and 
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1 observations would imply a convective turnover rate between 0.1- 0.2 day' 1 . However, if 

2 we assume the DC-8 had a positive bias toward sampling fresh convection in accordance 

3 with Fuelberg et al., our observed frequency distributions are consistent with a convective 

4 turnover rate closer to 0.1 day" 1 ( 45 , 58 ). 

5 

6 For comparison, the model detrainment cloud mass flux between 400 and 200 hPa 

7 (approximately 7.2 - 11.8 km based on standard atmosphere) for the INTEX-NA 

8 sampling domain (80° - 100°W for 30° - 35°N and 70° - 100° W for 35° - 50° N) between 

9 July 1 st and August 15 th 2004 was 0.0085 kg m' 2 sec' 1 . This corresponds to a turnover 

10 rate of 0.37 day" 1 (using a column mass of 1.9xl0 3 kg m' 2 between 7.2 - 11.8 km). As 

1 1 expected, the observed and modeled detrainment rates for the extratropics are slower than 

12 in the tropics. For example, Folkins and Martin determine a maximum convective 

13 detrainment rate of 0.4 day' 1 at 12.5 km from calculation of the clear-sky radiative and 

14 precipitation induced evaporative cooling rates using a one-dimensional model 

1 5 constrained by observations of temperature and water vapor ( 59 ). 

16 

17 Conclusions 

18 We present a statistical representation of the aggregate effects of convection on the 

19 chemistry and dynamics of the upper troposphere using in situ measurements taken 

20 aboard the NASA DC-8 during the summer of 2004 over the Eastern U.S. and Canada. 

21 These observations provide a new and unique constraint on: i.) the extent to which 

22 convection perturbs the continental UT during summer, ii.) the fraction of boundary layer 

23 air present in convective outflow, and iii.) the convective overturn rate of the upper 
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1 troposphere. In addition, the chemical clock described here defines a coordinate that can 

2 be used to assess the chemistry occurring down- wind of convective injection. These 

3 direct measures of atmospheric rates present a new opportunity for quantitative tests of 

4 model representations of processes governing UT ozone, convection, lightning and their 

5 impact on climate. 
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Upper Troposphere 


High NO x /HN 0 3 ratio 

0 5 10 IS 20 

Days After Convective Injection 

Addition of Lightning NO x 
* * Wet Scavenging of HN0 3 




Boundary Layer 


Figure 1: In moist convection, air from near the Earth’s Surface is rapidly transported 
upwards and detrained into the Upper Troposphere. In this process, Nitric Acid (highly 
soluble) is efficiently scavenged, while NO x (insoluble) remains. NO x is dramatically 
elevated by concurrent lightning NO production, resulting in high NO x to HNO 3 ratios in 
the convective outflow region. Following detrainment into the UT, NO x is converted to 
HNO 3 by OH during the day and via NO 3 /N 2 O 5 at night. The chemical evolution of the 
NO x /HN 03 ratio provides a unique indicator of the time a sampled air-mass has been in 
the UT following convection. 
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Figure 2: The observed deviation of the NO x to HNO3 ratio from steady-state as a 

5 function of altitude in the UT. The mean values within 500 m vertical bins are shown 

6 with circles (o). The steady-state NO x to HNO3 ratio was calculated from measured NO x , 

7 OH and Jhno3 and includes the rain out parameterization of Giorgi and Chameides 

8 (1985). 
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Figure 3: top panel Time series of measurements taken in the vicinity of recent 
convective activity on 1 1 August 2004 between 5 and 9 km. Panel A suggests the 
sampling of a series of fresh convective plumes, indicated by a sharp increase in the 
NO x /HN 03 ratio. Panels B & C depict coincident enhancements in SO 2 and UCN (3nm > 
D p > lOnm) and coincident sharp drops in CO 2 , indicative of the convective lofting of 
boundary layer air depleted in CO 2 . The derived time since the sampled air-mass had 
been influenced by convection is shown in Panel D. bottom panel NLDN lightning hits 
(small dots) on the 10 th and 1 1 th of August. The color-code represents the time of the hit 
(hours) prior to aircraft sampling. The DC-8 sampling location corresponding to 
measurements shown in Figure 1 is located on the Maine - New Brunswick border [46°N 
67°W]. The two day back trajectory [•] (initialized at the point of the second convective 
plume shown in Panel A) is also color-coded by time prior to DC-8 sampling (dots with 
black edges). 
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Figure 4: Observations of the fraction of ultra-fine condensation nuclei [number density 
of aerosol (3-10 nm) / total aerosol number density] (Panel A), NO x (Panel B), HNO 3 
(Panel C) and NO y (Panel D) as a function of modeled time since convective influence. 
The mean (-□-) and median (-0-) of the observations, within 8 hour bins, is shown along 
with interquartile range (shaded region). Results from the time-dependent box model, 
initialized at 10 km and 12PM, are shown with dashed lines for the gas phase species (B- 
D). 
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1 Days S trice Convective Influence 

2 

3 Figure 5: Observations of CO (Panel A) and Ozone (Panel B) as a function of modeled 

4 time since cloud processing in the UT. The mean (-□-) and median (-o-) of the 

5 observations, within 8 hour bins, is shown along with the interquartile range (shaded 

6 region). Results from the time-dependent box model, initialized at 10 km and 12PM, are 

7 shown with dashed lines. 
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Figure 6: left panel Normalized frequency distribution in the time since convective 
influence, as calculated from observations of the NO x to HNO3 ratio made during the 
summer of 2004. Calculations are separated into 1 km altitude bins (ranging from 7.5- 
11.5 km). The fraction of air that had been influenced by convection within the past two 
days (f < 2 days) is included in the figure legend, right panel Comparison of observed 
frequency distribution (7.5-11.5 km) with various modeled representations of the 
convective turnover rate. 
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1 Supplemental Online Information 

2 1. INTEX-NA Experiment Description and Instrument Descriptions 

3 The Intercontinental Chemical Transport Experiment - North America (INTEX-NA) 

4 took place between 1 July and 14 August 2004. Research flights were conducted out of 

5 Dryden Flight Research Center (Edwards AFB, CA), Mid-America Airfield (Mascoutah, 

6 IL); and PEASE International Trade -Port (Portsmouth, NH). Figure SI depicts the 

7 vertical and horizontal extent of research flights conducted aboard the NASA DC-8 

8 during INTEX-NA (7). DC-8 flight tracks are shown in the left panel of Figure 1 and the 

9 number of samples (10 second averaging time) in 1km vertical bins are shown in the right 

10 panel. In situ observations relevant to this study include; NO2, HNO3, OH, O3, CO, CO2, 

11 SO2 and Ultra-fine Condensation Nuclei (UCN). Table SI describes the detection 

12 threshold, uncertainty and time response for each measurement used in this analysis. 

13 

14 NO2 LIF Instrument The NO2 instrument flown aboard the DC-8 was described in 

15 detail by Thornton et al. (2), with specifics of the jet-expansion described by Cleary et al. 

16 (3). Briefly, NO2 fluorescence is detected at 1Hz following excitation of a specific jet- 

17 cooled rovibronic transition in NO2 at 585 nm. Red-shifted fluorescence is imaged at 90° 

1 8 onto an air cooled photomultiplier tube that is both optically and temporally filtered to 

19 remove laser scatter. NO2 mixing ratios are calculated directly from fluorescence counts 

20 following calibration to NO2 gas standards and measurements of the instrument zero from 

21 compressed air mixtures containing zero NO2. Selectivity to NO2 is demonstrated by 

22 tuning on and off of a specific NO2 resonance, where the difference in observed 

23 fluorescence at the two different frequencies is attributed solely to NO2. We calculate a 
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1 NO 2 detection threshold of 8 pptv in 10 seconds at the surface and 25 pptv in 10 seconds 

2 at the aircraft ceiling (12.5 km). 

3 

4 2. 0-D Time Dependent Model 

5 The chemical evolution of convective outflow was modeled using a 0-D time dependent 

6 model. The model was initialized with chemical conditions, altitudes and detrainment 

7 times consistent with observations of fresh convection made during INTEX-NA. As time 

8 propagates in the model, we calculate the production and loss of O 3 , CO, NO, NO 2 , NO 3 , 

9 N 2 O 5 , PAN, HO 2 NO 2 , HN0 3 , OH, H0 2 , R0 2 , H 2 0 2 , CH 3 OOH, H 2 CO and C1-C6 

10 Hydrocarbons for 20 days following cloud detrainment. The conversion of NO x to HNO 3 

11 in the outflow region is used as an indicator of time since convection. Figure S2 depicts 

12 the results of a single run initialized at 10km with a noon detrainment time. Initial 

13 conditions correspond to [NO x ]i = 800 pptv, [ 03 ]; = 65 ppbv and [CO]; =105 ppbv. Rapid 

14 conversion of NO x to HNO 3 is observed during the first few days as the system 

15 approaches steady-state. In this analysis we assume: i.) HNO 3 is scavenged with unit 

16 efficiency in deep convection, ii.) YN 205 = 0.01 and iii.) HNO 3 is not scavenged by 

17 aerosols (or rain) following injection into the UT. All kinetic rates used in this analysis 

18 were taken from the NASA JPL Chemical Kinetics and Photochemical Data for Use in 

19 Atmospheric Studies, Evaluation Number 14 ( 4 ). 

20 

21 2.1 Treatment of OH and HO 2 

22 The calculated time since convective detrainment is directly coupled to the HO x budget 

23 through the daytime NO x sink to HNO 3 via reaction with OH. As in other model 
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1 descriptions of the UT during INTEX-NA ( 5 , 6 ), our unconstrained model over-estimates 

2 OH by nearly a factor of two in the UT and under-estimates HO 2 by a similar amount. 

3 Due to the direct dependence of our timing indicator on HO x , we constrain the mixing 

4 ratios of OH and HO 2 to the observed values as a function of NO x and altitude. Figure 

5 S3 depicts the modeled mixing ratios of OH and HO 2 (lines), constrained to the 

6 observations (dots), as a function of NO x and SZA at 10km. The observed OH is a strong 

7 function of NO x , while observations of HO 2 remain insensitive to NO x . Constraints for 

8 OH and HO 2 were derived independently for each 1km altitude bin. Constraining OH 

9 and HO 2 to the observations increases the time required for the N 0 x -HN 03 system to 

10 reach steady-state (by slowing the rate of OH + NO 2 ) and enhances the modeled O 3 

1 1 production in the outflow region (by speeding up the rate of HO 2 +NO). 

12 

13 2.2 Calculation of Time since Convection 

14 The time since a sampled air-mass had been cloud processed is calculated by applying the 

15 mapping of time to NO x / HNO 3 derived in the box model to the observed NO x to HNO 3 

16 ratio. Figure S4 depicts the best-fit relation between the modeled NO x to HNO 3 ratio and 

17 time since cloud processing at 10km. This function is calculated at 1km increments from 

18 6-1 2km and applied to the measured NO x to HNO 3 ratio. 

19 

20 2.3 Model Assumptions and Uncertainty 

21 In order to access the uncertainty in the calculated time, we ran the time-dependent model 

22 under various different conditions encountered during INTEX-NA (e.g. [NO x ], (0. 2-3.0 

23 ppbv), [CE]; (40-100 ppbv), [CO]; (80-150 ppbv), detrainment time (noon, 4PM, 
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1 midnight), altitude (6- 12km) and time of year (June-September). As illustrated in figure 

2 S4, the NO x to HNO3 ratio has good resolution (i.e. large rate of change per unit time) 

3 during the first five days following convection. Beyond five days small changes in 

4 N0 x /HN 03 correspond to larger changes in the derived time. From the variance in the 

5 calculated time of individual model runs, we estimate the uncertainty in our modeled time 

6 to be ±6 hours at 1 day, ±12 hours at 2 days and ±1 day at 4 days. In addition, the 

7 INTEX-NA sampling domain did not permit frequent measurement of aged convection 

8 (>5 days). For these reasons we limit our analysis to the first five days following 

9 convection. 

10 

1 1 2.4 Treatment of Mixing 

12 The mixing rate was determined by iterating the model until we had closure between the 

13 observed and modeled time evolution of a suite of long-lived species (e.g. CO, CH4, 

14 CH3OH and others). The modeled time rate of change of species X is determined as: 

15 ^ =P(X) - L(x)-k ma (x] - [x } taclgro ,„ ) 

16 where P(X) represents the chemical production of species X, L(X) represents the 

17 chemical loss of species X and kdiiution is the mixing rate of the convective plume with 

18 background UT. We find this mixing term to be on average 0.05 ± 0.02 day' 1 . That is, 

19 after 5 days, the plume still has 75% of its original contents. While individual convective 

20 plumes may mix faster (or slower) that this, the aggregate mixing rate of all the sampled 

21 plumes can be described by this rate. This rate supports the conclusion that over the 

22 course of 5 days, the convective plume remains relatively isolated from the background 
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1 UT. Background mixing ratios used in the dilution calculation were taken as the mean 

2 observed UT mixing ratios outside of fresh convective plumes. 

3 

4 Due to subsidence of convectively lofted air parcels following injection, our calculated 

5 time represents a lower bound for age as the chemical clock speeds up (due to NO x 

6 repartitioning) as the parcel descends in altitude. However, this is a relatively small 

7 effect as calculated subsidence rates are approximately 35 hPa day' 1 (7). 

8 

9 3.0 Comparison of Chemical and Meteorological Convective Influence Calculations 

10 The results presented here provide a chemical constraint on the rate at which the UT over 

11 the continental US is influenced by convection during summer. In addition to the 

12 meteorological analysis of Fuelberg et al., presented in this manuscript, Thompson et al. 

13 assessed the effects of convection on the O3 budget in the UT, over eastern North 

14 America, during the summer of 2004. Using results from the INTEX Ozone Sonde 

15 Network Study (IONS), the authors conclude that 10-15% (lower-limit) of the below 

16 tropopause O3 can be attributed to the interaction of regional pollution with convection 

17 and lightning (5). 

18 

19 The INTEX-NA sampling period (June-August) and region (Eastern North America) is 

20 characterized by intense lightning activity. Cooper et al. calculated that 13% of the 

21 global annual lightning NO x emissions occurred between 108° W - 18 0 E and 18 0 -72° N 

22 the between June 21 and August 15, 2004 (9). In terms of lightning flash counts, 

23 Hudman et al., conclude that 2004 was typical (within 20% of the mean) of the past 5 
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1 years ( 6 ). Persistent frontal passages prevented the formation of stagnant high pressure 

2 systems, typically observed during the summer over the northeastern United States ( 10 ). 

3 These frequent passages led to both record low temperatures and number of O 3 

4 exceedances in the northeast ( 11 ). In contrast UT/LS O 3 , as observed from the IONS 

5 network, was comparable to the climatology ( 11 ). 

6 

7 4.0 Measurement Uncertainties 

8 In this analysis we calculate NO x from observations of NO 2 , O 3 , HO 2 and photolysis rates 

9 measured directly on the DC- 8 . NO was measured directly on the DC -8 via a 

10 commercial grade chemiluminescence detector. The sensitivity of the 

1 1 chemiluminescence instrument (detection threshold >50 pptv) and long integration time 

12 (1 minute) prohibited its use in these calculations. NO was calculated from steady-state 

13 using the following equation: 

[ at /}] = J _I2i \- N ° 2 ^ 

14 [i V uy J Steady-State 7 X ~\ 1 \ r 1 

k No + o 3 [0 3 \+k NO+H02 [H0 2 \ 

15 We calculate the accuracy of the derived NO to be better than ± 30% based on the 

16 propagation of the individual errors used on the calculation. Observations of Nitrogen 

17 Dioxide, Ozone and JNO 2 made during the INTEX-NA campaign were compared 

18 directly with measurements made aboard the NOAA WP-3D during a series of in-flight 

19 comparisons. During these experiments the principal individual components (NO2, O3 

20 and JNO 2 ) showed agreement to within their stated instrumental uncertainty. 

21 


6 



Bertram et al. Supplement 29 October 2006 Revision #2 


1 In Figure S6 we compare the measured NO x (using the chemiluminescence and the Laser 

2 Induced Fluorescence measurements) and NO x calculated from steady-state for the entire 

3 INTEX-NA campaign. The steady-state NO x agrees with the measured NO x to within the 

4 calculated uncertainty when averaged to 1 -minute and divided into 1km vertical bins. 

5 The observed upper tropospheric NO x concentrations during INTEX-NA are on average 

6 much higher than previous reported on intensive aircraft field campaigns over North 

7 America. During the SUCCESS campaign, Jaegle et al. report mean NO concentrations 

8 of 0.030 ± 0.022 ppbv and 0.061 ± 0.045 ppbv, for the altitude ranges of 8-10 km and 10- 

9 12 km, respectively (12). However, it must be noted that the scientific objective of 

10 SUCCESS was the sampling of aircraft exhaust and contrails, thus the values reported in 

11 Jaegle et al., were filtered to exclude fresh aircraft exhaust (CO 2 > 368 ppmv and NO > 

12 300 pptv). Crawford et al., report mean NO concentrations of 0.1 ppbv for observations 

13 made between 6-12 km during the SONEX campaign during the fall of 1997 over the 

14 North Eastern United States and the North Atlantic (13). Neither of these studies 

15 provides a direct comparison to the INTEX-NA data-set as SUCCESS was conducted in 

16 the spring and SONEX in the fall, while peak lightning and convection occurs over the 

17 continental North America during summer. The best comparison comes from NO x 

18 measured aboard a co mm ercial passenger aircraft during the NOXAR program between 

19 1995 and 1997. Brunner et al. report a mean UT NO x concentration of 0.4 ppbv during 

20 June, July and August over North America (14). As seen in Figure S6, our observations 

2 1 during the summer of 2004 are consistent with this. 

22 
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1 We use the Caltech CIMS HNO3 due to its fast time response (5 seconds as compared to 

2 105 seconds for the UNH Mist Chamber Technique) and the UNH MC results when the 

3 fast HNO3 was unavailable. To account for the systematic bias between the two 

4 observations in the UT ([HN03 ]unh = 0.6 x [HN03]c a itech), we scale both the CIMS and 

5 MC observations to split the difference between the two measurements (i.e. we increase 

6 [HN0 3 ]unh by 20% and decrease [HN 03 ]cit by 20%). 

7 

8 Due to the observed systematic bias, the choice of which HNO3 measurement to use in 

9 the analysis has the potential to complicate our conclusions. To address these effects we 

10 have conducted the analysis using Nitric Acid as measured from: i.) the Mist Chamber - 

1 1 Ion Chromatography Instrument, ii.) the Chemical Ionization Mass Spectrometer and iii.) 

12 the scaled difference between the two techniques. The results are compared in Figures 

13 S7 and S8. Figure S7 shows the normalized frequency distribution in the time since 

14 convective influence, as calculated from observations of the NO x to HNO3 ratio made 

15 during the summer of 2004. Figure S8 depicts the fraction of air that had been influenced 

16 by convection within the past two days (f < 2 days) as a function of altitude. Calculations 

17 derived from the Mist Chamber - Ion Chromatography Instrument result in a higher 

18 fraction of convectively influenced air when compared with calculations made using 

19 measurements from the CIMS instrument. When compared with the results shown in 

20 Figure 6B of the manuscript, on which our conclusions regarding the convective overturn 

21 rate are drawn from, calculations using either the CIMS, MC-IC or the scaled difference 

22 result in a convective overturn rate between 0.1 and 0.2 day' 1 . 
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1 2. Supplemental Figures 



140°W 120°W 100°W 80°W 60°W 40°W o iu m in a ■ u> m 

Number of Samples. (Using 1 0-see Dala} 


3 Figure SI: left panel INTEX-NA flight tracks made between 1 July 2004 and 14 August 

4 2004 aboard the NASA DC-8, right panel Number of samples (using 10-sec averaged 

5 data) within 1km altitude bins between 0-12 km during the entire campaign. 
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Days Since Convection 


3 Figure S2: Time-dependent model illustrating the conversion of NO x to Nitric Acid in 

4 the days subsequent to a cloud processing event occurring at 10km. The above model 

5 was initialized at 12PM local time at 30°N in August using [NO x ]; = 800 pptv, [CO]i = 

6 105 ppbv at [ 03 ]; = 65 ppbv. 
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4 Figure S3: Model representation of OH (left panel) and HO 2 (right panel) as a function 

5 of SZA and [NO x ]. Model results (solid lines) are shown on top of the in situ 

6 observations (dots). The model was initialized at noon at 10km with [NO x ]i = 800 pptv, 

7 [CO]; = 105 ppbv at [ 03 ]; = 65 ppbv. Observations shown were taken aboard the DC -8 

8 between 9 and 1 1 km. 


11 


Men 


Bertram et al. Supplement 29 October 2006 Revision #2 


1 



[NOJ / [HN0 3 ] 


5 Figure S4: Observed NO x to HNO3 ratios are converted to a time since last convective 

6 influence using the best fit equation relating the N 0 X /HN 03 ratio to time as calculated 

7 using the time-dependent model in 1km altitude bins from 7.5-11.5 km. The above 

8 equation is valid for pressure altitudes between 9.5 and 10.5 km. 
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Days Since Convective Influence Influence (or> Sdays) 

1 

2 

3 Figure S5: Comparison of chemical (grey bars) and meteorological constraints (-o-, -□-) 

4 on convective influence during INTEX-NA. Convective influence on air sampled by the 

5 DC-8 is shown with blue circles (-o-), while convective influence on the entire fNTEX- 

6 NA domain is shown with red squares (-□-). 

7 

8 
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2 

3 Figure S6: Comparison of observed (grey lines) and steady-state (black lines) NO x for 

4 the entire INTEX-NA field campaign. The shaded region represents the interquartile 

5 range of the calculated NO x . 

6 
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Days Since Convective Influence Influence (or > 5 days) 

2 

3 

4 Figure S7: Normalized frequency distribution in the time since convective influence, as 

5 calculated from observations of the NO x to HNO3 ratio made during the summer of 2004. 

6 Calculations were conducted using Nitric Acid as measured from: i.) the Mist Chamber - 

7 Ion Chromatography Instrument (-□-), the Chemical Ionization Mass Spectrometer (-0-) 

8 and the scaled difference (-■-). 
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2 

3 Figure S8: Fraction of air that had been influenced by convection within the past two 

4 days (f < 2 days) as a function of altitude. Calculations were conducted using Nitric Acid as 

5 measured from: i.) the Mist Chamber - Ion Chromatography Instrument (-□-), the 

6 Chemical Ionization Mass Spectrometer (-0-) and the scaled difference (-■-). 

7 

8 
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1 


Species 

Measurement 

Technique 

Detection 

Threshold 

Accuracy 

Time 

Response 

Reference 

no 2 

LIF 1 

8 pptv / 10 sec 

± 10% la 

1 Hz 

(2,3) 

hno 3 

CIMS 3 

10 pptv / 0.5s 

± 30% 

0.5s 
sample 
every 5 s 

( 15 ) 


Mist Chamber - IC 4 

5 pptv / 105 sec 


105 sec 

( 16 ) 

OH 

LIF 1 

0.01 pptv 

± 32% 
2a - 1 min 

20 sec 

( 17 ) 

o 3 

Chemiluminescence 

Precision = 
±0.8pptv,±l% 

of reading 

± 2 ppbv, 
±3% 

1 Hz 

( 18 ) 

CO 

IR- Absorption 

Precision = 
±lppbv, ±1.5% 
of reading 

±1.4 ppbv, 
± 2.6% 2a 

1 Hz 

( 19 ) 

co 2 

IR- Absorption 

Precision < 
0.07 ppmv 

± 0.25 ppmv 

1 Hz 

( 20 ) 

UCN 

TSI CN counter 5 

± 10% 

± 10% 

1 Hz 

( 21 ) 

Jn02 

Actnic Flux 
Spectroradiometer 

4.1 x 10' 7 
Precision = 
±4.5% 

±8% 

1 Hz 

( 22 ) 


2 'LIF - Laser Induced Fluorescence 

3 N0 2 detection threshold is 8 pptv / 10 sec at 760 Torr (ground) and 20 pptv / 10 sec at 200 Torr (10 km) 

4 3 CIMS - Chemical Ionization Mass Spectrometry 

5 4 IC - Ion Chromatography 

6 5 UCN (Ultra-fine Condensation Nuclei) was obtained by the difference of the UCN (D p >3nm, TSI 3025) 

7 and CN (D P >10nm, TSI 3010) Condensation Nuclei (CN) instruments. 

8 

9 Table SI: Detection thresholds, measurement uncertainty and time response of the in 

10 situ measurements used in this study. (Note: ppmv = pmol mol' 1 , ppbv = nmol mol' 1 and 

1 1 pptv = pmol mol' 1 ) 
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